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ABSTRACT 

The planetary nebula TS01 (also called PNG 135.9+55.9 or SBS 1150+599A), with its record-holding low oxygen abundance and 
its double degenerate close binary core (period 3.9 h), is an exceptional object located in the Galactic halo. 

We have secured observational data in a complete wavelength range in order to pin down the abundances of half a dozen elements in 
the nebula. The abundances are obtained via detailed photoionization modelling taking into account all the observational constraints 
(including geometry and aperture effects) using the pseudo-3D photoionization code Cloudy _3D. The spectral energy distribution of 
the ionizing radiation is taken from appropriate model atmospheres. Incidentally, from the new observational constraints, we find that 
both stellar components contribute to the ionization: the "cool" one provides the bulk of hydrogen ionization, and the "hot" one is 
responsible for the presence of the most highly charged ions, which explains why previous attempts to model the nebula experienced 
difficulties. 

The nebular abundances of C, N, O, and Ne are found to be respectively, 1/3.5, 1/4.2, 1/70, and 1/11 of the Solar value, with uncer- 
tainties of a factor 2. Thus the extreme O deficiency of this object is confirmed. The abundances of S and Ar are less than 1/30 of 
Solar. The abundance of He relative to H is 0.089±0.009. 

Standard models of stellar evolution and nucleosynthesis cannot explain the abundance pattern observed in the nebula. To obtain an 
extreme oxygen deficiency in a star whose progenitor has an initial mass of about 1 M e requires an additional mixing process, which 
can be induced by stellar rotation and/or by the presence of the close companion. We have computed a stellar model with initial mass 
of 1 M , appropriate metallicity, and initial rotation of 100 km s -1 , and find that rotation greatly improves the agreement between the 
predicted and observed abundances. 

Key words. ISM: planetary nebulae: individual — ISM: abundances — Stars: AGB and post-AGB - Stars: binaries: general — 
Physical data and processes: Nuclear reactions, nucleosynthesis, abundances 
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1. Introduction 

SBS 1150+599A was discovered in the second Byurakan 
Sky Survey and firs t classified as a c atacly smic variable 
(IStepanian et al.l Il999h . iTovmassian et al.l (l200lh discussed in 
detail the nature of the object and came to the conclusion 
that, in fact, it is a planetary nebula (PN). The object was 
renamed PN G 135.9+55.9, following the nomenclature for 
Galactic PNe from the S trasbourg-E S O catalogue of Galactic 
Planetary Nebulae ( I Acker et al.lll992l) . For the sake of brevity, 
we will refer to it as TS01 in the rest of the paper. This PN 
is special in at least three important aspects. First of all, its 
oxygen abundance is very low, significantly lower than in any 
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other PN known up to now ( Tovmassian et al. 2001; R icher et"aT1 
120021: IJacobv et al.l l2Q02t IPequignot & Tsamisll2005l) . Second, 
its nucleus is a spectroscopic bina ry, with a period of only a 
few hours (iTovmassian et al.ll2004h . Third, it appears, from es- 
timates of the nature and masses of the two stellar compo- 
nents, that TS01 could turn into a double degenerate Type la 
Supernova dTovmassian et al.ll2004l) . Each of these aspects, even 
taken alone, makes TS 01 an exceptional object. 

In this paper, we reexamine the chemical composition of 
TS01. Briefly, the story of the determination of the chemi- 
cal composition of this object is the following. Tovmassian 
et al. (2001) had optical spectra of TS01 in the range 3900- 
7000 A obtained with 2 m class telescopes which showed no 
lines from heavy elements except a very weak [O m] /15007, 
with an intensity a few percent of Hfi. A coarse photoion- 
ization analysis suggested an oxygen abundance smaller than 
1/100 of Solar. Note that standard empirical methods for abun- 
dance determinations in PNe cannot be used for TS01, since 
the electron temperature cannot be determined directly from 
observations. To go further in the abundance determination of 
TS01 required an estimate of the effective temperature of the 
central star. One way is to obtain a good blue spectrum of 
the PN, and use the [Ne v l ,13426/ [Ne ml ,13 869 ratio (or a 



limit on it) as a constraint. iRicher et al 



J2002I) at th e Cana da- 



France-Hawaii Telescope (CFHT) and iJacobv et al.l (120021) at 
the Multiple Mirror Telescope (M MT) secured deep b lue spec- 
tra in order to detect these lines. IJacobv et al. (12002 ) d etected 
the [Ne v] ,13426 line at a level of 0.8 H/J. IRicher et al.l d2002l) 
found only an upper lim it of 0. 1 H/3! Concerning the [Ne m] 
/13869 line, IJacobv et al.l d2002l) measur ed an intensity about 
10 times larger than IRicher et al.l d2002l) . The two papers ap- 
peared within a few days of each other on astro-ph, revealing 
this big conflict in the observations. The two groups conducted 
independent photoionization analyses, and both concluded that 
the O/H ratio is less than 1/100 of Solar (the main reason for 
their similar result for the oxygen abundance was the simi- 
lar [Ne v] ^3426/[Ne ml ,13869 ratio used by both studies). 
IPequignot & Tsamisl (120051) merged and discussed the two ob- 
servational data sets and conducted their own photoionization 
analysis. They concluded that the O/H ratio of TS01 lies be- 
tween 1/30 - 1/15 of Solar (still holding the record for the most 
oxygen poor planetar y nebula but much hi g her th an previously 
published). However. IPequignot & Tsa mis (2005) neglected to 
consider observations of TS01 made with the Hubble Space 
Telescope (HST) and the Far Ultraviolet Spectroscopic Explorer 
(FUSE). As a result, some of their "predicted" line intensities 
are in conflict with what is actually observed in the UV. HST ob- 
servations were obt ained in 2003 , and p resented in a short, pre- 
liminary version by IJacobv et al.l (120061) . Those authors quoted 
an oxygen abundance of 1/30 - 1/40 of Solar, and carbon and 
nitrogen abundances roughly 1/10 of Solar. 

Before embarking on a new determination of abundances, 
we have chosen to gather the best possible observations at all 
wavelength ranges. These data provide many more constraints 
than were available in any previous study. In order to make the 
best use of the large amount of data obtained with different tele- 
scopes, we use a pseudo-3D photoionization code, Cloudy _3D, 
which is able to account for the nebular geometry as we see 
it now, and with which we can properly take int o account the 
apertu re effects. This cod e is based on CL OUDY dFerland et alj 
1 1998b and was written bv lMorissetl d2006l) . 

The paper is organized as follows. Section 2 presents the new 
observational material: several optical spectra, HST imaging and 
spectroscopy, infrared spectroscopy with the Spitzer Telescope, 



and mentions our X-ray observations with XMM. Section 3 sum- 
marizes other data that we used as constraints for the photoion- 
ization modelling. Section 4 describes our modelling strategy, 
and presents our "reference model". Section 5 evaluates the error 
bars on the derived elemental abundances, taking into account 
observational uncertainties in emission-line fluxes, uncertainties 
in model input parameters as well as uncertainties arising from 
an imperfect description of the physical processes included in 
the models. In Section 6, we compare the chemical composition 
of TS 01 with that of other PNe in the Galactic halo and discuss 
it in terms of stellar nucleosynthesis in the Asymptotic Giant 
Branch (AGB) phase. Finally, Section 7 summarizes our main 
findings. 



2. New observational data on emission lines 

We present the observational data that we secured on TS 01 and 
its stellar core since the work presented in Tovmassian et al. 
(2004). Some of those data were already briefly reported in con- 
ference proceedings, but here we describe the acquisition and 
reduction processes in more detail. Note that the observations 
and reductions were done by different people and at different 
epochs, when our knowledge on the object was not the same. 
This explains the differences in the tactics employed to reduce 
the data, estimate the line fluxes and correct for reddening. We 
did not try to fully homogenize the data reduction process, since 
we felt it unnecessary for our purposes. 

The log and characteristics of each set of observations are 
given in Table Q] 



Table 1. Log and characteristics of the spectroscopic observa- 
tions 



telescope 


date 


,1 range 


resolution 


aperture 


FUSE 


30 Jan 02 


900-1200 A 




30x30" 


HST STIS 


4 May 03 


1170-1700 A 


1.20 A 


0.5" 


CFHT MOS 


1 May 03 


3400-5300 A 


3.0-3.5A 


1" 


CFHT MOS 


4 Mar 01 


3400-8000 A 


23A 


5" 


Kitt Peak 


1 Jan 03 


3600-7500 A 


7A 


2" 


SDSS 




3819-9196 A 


2-4A 


diameter 3 " 


Spitzer IRS SH 


22 Apr 06 


9.9-19.6/mi 


600 


4.7 x 11.3" 


Spitzer IRS LH 


22 Apr 06 


18.7-37.2pm 


600 


11.1 x 22.3" 



2.1. Imaging 

We (M.P.) retrieved from the HST archives and analyzed the data 
corresponding to the proposal ID 9466. The observations were 
performed on May 5, 2003. Two types of data are available: di- 
rect imaging and spectroscopy. 

Direct imaging was obtained with the Advanced Camera for 
Surveys (ACS). The High Resolution Channel with a field of 
view of 26"x29"and a plate-scale of 0.027"per pixel, with filters 
around He (central wavelength 6581.97 + 162.8 A) and [Ne v] 
,13426 (central wavelength 3432.8 + 42.66 A) were used. 

Figure Q] shows the Ha image obtained by averaging 
the 4 calibrated frames: j8do01021, j8do01022, j8do01023, 
j8do01024 (870 s exposure time each; 58 min in total), after 
aligning them with respect to j8do01021. The image is roughly 
elliptical in shape with two brighter, symmetrically-placed lobes 
at a position angle of about 103° that extend the full major axis 
of the ellipse. The nebula is not perfectly symmetric, with the 
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Fig. 1. HST-ACS Ha image of TS 01, with a logarithmic grey Fig. 3. HST-ACS image of TS 01 in the [Ne v] ,13426 line, on a 
scale. A grid with the position scale (R.A. and Dec) is traced, logarithmic grey scale. Same orientation and scale as Fig.Q] 
North is up. 




Fig. 2. HST-ACS Ha image of TS 1 . The locations of the differ- 
ent spectroscopic apertures are indicated: Kitt Peak (continuous 
lines), HST STIS (dotted lines), SDSS (circle). See text for the 
CFHT 2003 observations. 



outermost southern part much fainter. The size of the nebular 
image is about 5". 

Figure [2] shows the same image as Fig. Q] with the differ- 
ent observing apertures indicated: continuous lines for Kitt Peak, 
dotted lines for HST STIS, circle for SDSS. For the CFHT 2003 
observations, the slit was rotated before each of the seven ex- 
posures (see Sect. 12.2.11 ), so as to remain as close as possible 
to the parallactic angle. The position of the slit is not indicated 
in the figure, for the sake of clarity, but was taken into account 
correctly when comparing the predicted line intensities with the 
observed ones (see Sect. 14.21 ). 

Figure [3] shows an average of the [Ne v]/l3432 calibrated 
images j 8do0 1 1 1 , j8do0 1 1 2, j8do0 1 1 3 and j 8do0 1 14 (3000 
sec of exposure time each; 200 minutes in total). This image re- 
veals a very faint, roughly spherical extended nebulosity and an 
important emission in the center, probably caused by the stellar 
core. Some faint extensions are marginally detected in the direc- 
tions of the Ha lobes. 



2.2. Optical spectroscopy 
2.2.1. CFHT data 

CFHT 2003 TS01 was reobserved at the CFHT by M.R. and 
G.S. on 1 May 2003 using the MOS spectrograph and a 1" slit 
(lLeFevreetal.1119941) . The U900 grism was used, giving a spec- 
tral range of 3400-5300 A and a spectral resolution of 3-3.5 A 
(measured from arc lamp spectra). Seven 1800 s exposures were 
obtained. During each exposure, the slit was set to within 10° 
of the parallactic angle. Details of the reduction process of the 
individual exposures are given in Tovmassian et al. (2004). 

To obtain a high signal-to-noise spectrum of the nebular 
emission lines, it is necessary to account for the stellar and neb- 
ular continuum emission. These contributions were subtracted 
from the individual exposures before summing the individual 
spectra. First, the observed spectra were shifted in velocity so 
that the stellar absorption line was at rest. Next, the H/3 intensity 
was measured and used to scale a model of the nebular contin- 
uum emission. Since the Hp flux is affected by stellar absorption, 
the stellar absorption was assumed to have an equivalent width 
of 13 A, a value typical for the models used (see below). This 
scaled nebular continuum was then subtracted from the observed 
spectrum. Then, a model stellar atmosphere was scaled so as 
to match the observed continuum level and subtracted from the 
observed spectrum. This procedure leaves a pure nebular emis- 
sion line spectrum, supposing that the model nebular and stellar 
continua are representative of their real counterparts. It is un- 
likely that subtracting the continua introduces significant uncer- 
tainty into our final line intensities. Model stellar spectra with 
(T eff ,logg) pairs of (90 kK, 5.05), (120 kK, 5.35), and (150 kK, 
5.56) computed by R.N. (see Tovmassian et al. 2004) were sub- 
tracted from our observed spectra and the differences in the re- 
sulting line intensities were always smaller than the uncertainties 
in the fitfl 



1 This treatment was applied well before we had understood that the 
optical continuum was dominated by a star of 55 kK (see Sect. [33}. In 
view of the fact that the lines we use for the dignostics discussed in the 
present paper are hardly affected by this correction, we decided not to 
redo the subtraction using more adequate model stellar spectra. 
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Fig. 4. Flux-calibrated CFHT spectra of TS01. In blue is the 
CFHT 2003 spectrum with the nebular and stellar continua sub- 
tracted, in black the original CFHT 2001 spectrum, and in red 
the CFHT 2001 spectrum with the correct wavelength solution. 



Once the stellar and nebular continua were subtracted, the 
individual nebular spectra were shifted back to their original ve- 
locities and summed. We measured the nebular emission line 
strengths from this final spectrum. The line intensities were mea- 
sured using INTEN S, a locally-implemented software package 
dMcCall et al.lll985l) . This software simultaneously fits a sam- 
pled Gaussian function to the emission line(s) and a straight 
line to the continuum. It returns the line strengths, line wave- 
lengths, and uncertainties in these quantities. The line intensities 
presented in Table 2, together with their uncertainties, are those 
measured after subtracting the stellar spectrum for r e ff = 120kK 
and log g = 5.35. The listed intensities are not corrected for red- 
dening. In the case of no detection, two-sigma upper limits are 
given instead. 

We note that the [Ne v] A3426 line is present, and strong. 
Its intensity is of the same order as found in the spectrum 
of iJacoby et al.l (120021). and much higher than the upper limit 
given by Richer et al. d2002l) . The remaining lines have intensi- 



ties ro ughl y in agreement with those published by iRicher et al.l 
d2002l) and ljacobv et al.l d2002l) . ex cept for the [Ne m] l ine which 
appeared on the top of a bump in IJacoby et al.1 (|2002|) and was 
attributed a high intensity in that paper. 



CFH T 200 1 In v iew of the large discrepancy with the 
IRicher et al.l (I2002 ) data concerning the [Ne v] A3A26 line, we 
decided to reanalyze the spectrum of TS 01 we had obtained in 
March 2001 at the CFHT. First, though, we refer the interested 
reader to Richer et al. (2002) for a discussion of the details of 
these observations. The basic difficulty with these observations 
was that the arc lamp spectra were taken with the same 5"slit 
used for the object observations, which resulted in an arc spec- 
trum with severely-blended lines. After repeating the wavelength 
calibration more carefully, we found that our previous solution 
had stretched out the spectrum at the shortest wavelengths, lead- 
ing us to not recognize the [Ne v] A3426 line because of its low 
contrast with respect to the continuum and its erroneous wave- 
length. The analysis of this spectrum was considerably simpler 
given that the stellar features were not resolved. We simply fit the 
continuum shape and removed it (INTENS assumes that the con- 
tinuum is a straight line), then measured the line intensities with 
INTENS. The resulting line intensities are given in Table 2. The 



intensity we find now for [Ne v] A3426 is in agreement with the 
one given by IJacoby et al.l d2002l) . and compatible within 2 sig- 
mas with the one obtained with the CFHT 2003 data. Slit effects 
could perhaps explain the slight difference in intensity between 
CFHT 2001 and CFHT 2003. 



2.2.2. Kitt Peak data 

TS 01 was observed at the 4 m telescope of Kitt Peak National 
Observatory on 1 January 2003. The grating used was KPC-10A, 
and the slit 2"x300", with an orientation of P.A.=44.7°. Two ex- 
posures of 600 s were obtained. The data were reduced by L.J., 
employing the same procedure as for the SDSS spectrum, ex- 
plained below. 

2.2.3. SDSS and Kitt Peak data 

The spectrum of TS 01 appears in the data of the Sloan Digital 
Sky Survey SDSS (http://www.sdss.org) under the name 0953- 
5241 1-16Cq We present here its analysis as performed by L.J. 

We separated the nebular emission from the stellar spectrum 
and evaluated the reddening with as few free parameters as pos- 
sible. We assumed the stellar spectrum to be that of a single 
white dwarf (WD), hence neglecting the possible contribution of 
the companion; we considered three model WD spectra at tem- 
peratures of 90, 120 and 150 kK (the same as used for the CFHT 
2003 spectrum). As for the nebular continuum, we computed 
the free-free and free-bound emiss ivities of H + and He ++ with 
the chianti code dLandi et al.ll2006l) . assuming an abundance ra- 
tio He ++ /H + = 0.075 and an electronic temperature of 30 kK 
(Richer et al. 2002). We also retrieved the Hy nebu lar emissivity 
at this temperature from lStorev & Hummerld 19951) . 

First, we computed a model of the total (stellar+nebular) 
spectrum around the Hy line. In each of the WD spectral mod- 
els, the Hy line has a Voigt profile with a given equivalent width 
(EW) W wc j, Gaussian width <x wd and Lorentzian width a wd . As 
for the nebular emission, we computed the EW W ne b of the emis- 
sion line with respect to the nebular continuum. Furthermore, 
we assumed the real width of the nebular line to be much 
smaller than the instrumental one, so its observed Gaussian and 
Lorentzian widths, respectively <x inst and fli nst , are representative 
of the instrumental PSF Consequently, the observed widths of 
the stellar line are cr» = (cr^ d + cr? nst )^ 2 and a, = a wd + fli nst . We 
normalized the local stellar+nebular continuum with the fit of a 
slope on either side of the line. Finally, we let the central wave- 
lengths of the stellar and nebular line, respectively A* and A Be fo, 
be independent from each other. Using the data and assumptions 
gathered, we fitted a consistent model on the observed spectrum 
around the line. Calling V ne b(A) the profile of the nebular line, 
V*(A) that of the star (both being normalized to an EW of unity), 
and C* the stellar contribution to the flux, the model can be writ- 
ten as 

F A (A) = C* (1 -W wd V«(/l -A*; cr inst , a inst , <x wd ,a wd )) + 

+(1 - C) (1 + W neh V neb (A - A aeh ; cr inst , a inst )) (1) 



2 The analysis by L.J. was done soon after our discovery 
of the spectrum in SDSS d ata release 2. In data release 6 
jAdelman-McCarthv et al1l2008l) , SDSS spectra were recalibrated, re- 
sulting in an increase of about 30% of the fluxes of TS 01. Line ratios 
remained unchanged. Therefore, the analysis of L.J. remains valid. On 
the other hand, whenever we needed to consider the TS 01 continuum 
in this paper, we used the recalibrated spectrum. 
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Fig. 5. Fully processed SDSS spectrum (dereddened, free of con- Fig. 6. Fully processed Kitt Peak spectrum (dereddened, free of 
tinuum and stellar lines). continuum and stellar lines). 



For all three WD model atmospheres, we obtained very good fits 
with no visible systematic residuals. 

Given that the intrinsic Ha/HB nebular line ratio in TS01 
is not merely the case B recombination value, and given the 
additional problems with Ha (see Sect. I5.3.1I ), we cannot use 
this ratio to evaluate the extinction. Hence, we used the contin- 
uum to measure the latter. We first corrected the data for the 
small foreground ext inction {E(B - V) = 0.029) estimated by 
ISchlegel et all (119981) . Then, we removed most of the nebular 
or stellar lines from the observations making use of a median fil- 
ter. Finally, adopting the SMC extinction law dPrevot et al.l 19841 
lBouchetetai1ll985l) and comparing the filtered spectrum with 
the theoretical continuum, we evaluated the reddening and cor- 
rected the data for it. The choice of the extinction law was mo- 
tivated by the low metallicity of TS 01. The reddening amounts 
obtained (additionally to the foreground one) are E(B - V) = 
0.033, 0.044 and 0.050 for the 90, 120 and 150 kK WD models, 
respectively. 

The last processing of the data was the removal of the stel- 
lar and nebular continua, in particular to avoid the contamina- 
tion of the nebular lines by the underlying stellar features. We 
used the fit of the Hy line to shift the theoretical stellar spectrum 
and nebular continuum according to their evaluated radial veloc- 
ities, convolved them by the average instrumental PSF and sub- 
tracted them. Finally, we identified visually the detectable lines 
and measured their fluxes. The fully processed SDSS spectrum 
where the 90 kK WD model spectrum was removed is presented 
in Fig. [5] while Fig.[6]shows the result of the processing of the 
Kitt Peak spectrum. The choice of the WD model had a moder- 
ate impact on the evaluation of the line fluxes, of order of 2% 
for most of them and 2cr in the worst case. The intensities of the 
SDSS spectrum and the Kitt Peak spectrum are listed in Table|2] 

2.3. Ultraviolet spectroscopy 
2.3.1. HST STISdata 

The HST STIS spectroscopic data correspond to the same pro- 
posal (ID 9466) as the imaging data. 

A 52 " x 0.5 " slit was used. It was oriented along the bright 
jet-like emission of the nebula (PA. 103°), see Fig. [2] 

Far UV observations The MAMA detector combined with 
a G140L grating provided 2D spectra o8do03020, o8do03030, 




Fig. 7. 2D calibrated HST FUV spectrum of TS 01 



o8do03040, o8do03050 and o8do03060, with 4675 sees expo- 
sure time each. The spectra cover a wavelength range from 
1170 to 1700 A, with a resolving power of 1190 at the central 
wavelength 1425 A. The 2D spectra show a bright blue stellar 
continuum and a few faint and extended emission lines from 
the nebula. Calibrated 2D spectra were combined (after shift- 
ing because the spectroscopic observations were dithered) to 
produce a 389.6 min spectrum. The resulting 2D spectrum is 
shown in Fig. [7] The stellar spectrum shows good signal-to-noise 
and stellar and interstellar absorption are present. The analy- 
sis of the stellar spectr um is presented in a companion paper 
dTovmassian et all 12009). Regarding nebular lines, the follow- 
ing ones are detected: Nv^l240, Niv]/H488, Civ] ,11550 and 
He 11] A 1640. Selective absorption of resonance lines by the in- 
tervening interstellar medium is treated in Sect. 13.41 

From the combined 2D spectrum the nebular emission was 
extracted on both sides of the central star, with an extraction win- 
dow of 60 pixels, equivalent to 1.464". Fig. [8] is a combination 
of both nebular spectra. The line fluxes in each lobe, and the 
combined values with respect He n] A 1640 are listed in Table |3], 

Near UV observations The MAMA detector combined with 
a G230LL grating provided 2D spectra o8do02010, o8do02020, 
o8do02030, o8do02040 and o8do02050, covering a wavelength 
range from about 1600 to 3 150 A. The calibrated 2D spectra 
were combined (after aligning) to produce a spectrum with to- 
tal exposure time of 237.5 min. As for the FUV, the NUV stel- 
lar spectrum has good signal-to-noise and stellar and interstellar 
absorption can be seen. However, no nebular lines are detected. 
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Table 2. Intensities of optical lines, corrected for stellar absorp- 
tion, but not for reddening, with respect to H/3=100 





SDSS 


Kitt Peak 


CFHT 2003 


CFHT 2001 


[Ne V] 3426 






55.54 ± 5.19 


84.05 ± 13.07 


[O II] 3727 






< 0.60 


< 1.55 


H I 3735 






1.07 ± 0.31 




HI 3751 




4.94 ± 1.18 


1.70 ±0.32 




H I 3772 




4.36 ± 1.30 


2.30 ± 0.35 




H I 3798 




< 3.0 


2.64 ± 0.38 




H I 3836 


6.73 ± 0.53 


3.71 ± 0.82 


4.67 ± 0.37 




[Ne III] 3869 


< 1.2 


< 1.5 


0.77 ± 0.31 


1.67 ± 0.48 


H I 3889 


9.06 ± 0.43 


7.81 ± 0.87 


9.28 ± 0.45 


2.37 ± 0.47 


H I 3970 


13.29 ± 0.44 


13.24 ± 0.79 


12.70 ± 0.47 


5.85 ± 0.58 


He II 4027 


< 1.2 


< 1.5 


0.56 ± 0.21 




C III 4069 


2.1 :: 


< 1.5 


< 0.58 




HI 4102 


25.11 ±0.42 


25.49 ± 0.62 


25.16 ± 0.43 


19.79 ± 1.06 


He II 4201 


< 1.2 


< 1 


1.04 ±0.19 




C II 4267 


< 1.2 


< 1 


<0.32 




H 1 4340 


46.09 ± 0.38 


48.06 ± 0.44 


47.19 ± 0.50 


42.49 ± 0.80 


[O III] 4363 


< 1.2 


< 1 


< 0.41 




N III 4379 


< 1.2 


< 1 


< 0.09 




He 1 4471 


< 1.2 


< 1 


< 0.11 




He II 4543 


2.97 ± 0.33 


2.91 ± 0.35 


2.46 ± 0.22 




[Ar V] 4626 


< 1.2 


< 1 


< 0.38 




O IV 4632 


< 1.2 


< 1 


< 0.38 




C III 4650 


< 1.2 


< 1 


< 0.38 




C IV 4659 


< 1.2 


< 1 


< 0.38 




C IV 4659 


< 1.2 


< 1 


< 0.24 




He II 4686 


75.13 ±0.31 


79.02 ± 0.28 


77.50 ± 0.95 


79.06 ±1.11 


[ArlV] 4711 


< 1.2 


< 1 


< 0.15 




[Ne IV] 4715 


< 1.2 


< 1 


0.32 ± 0.08 




[Ne IV] 4725 


< 1.2 


< 1 


0.26 ± 0.08 




[Ar IV] 4740 


< 1.2 


< 1 


< 0.15 




H/3 4861 


100.00 ±0.29 


100.00 ± 0.26 


100.00 ± 0.78 


100.00 ± 1.59 


O V 4930 


< 1.2 


< 0.8 


< 0.28 




N V 4945 


< 1.2 


< 0.8 


< 0.28 




[O III] 4959 


< 1.2 


< 0.99 ± 0.29 


0.59 ± 0.13 


0.69 ± 0.69 


[O III] 5007 


1.82 ±0.28 


2.16 ± 0.26 


2.53 ±0.16 


1.81 ±0.79 


[FeVI] 5146 


< 10.0 


< 0.8 


<0.28 




Hell 5411 


5.8 ±0.22 


6.49 ± 0.28 




5.12 ± 0.31 


He I 5876 


< 0.8 


< 0.8 




< 0.19 


[Fe VII] 6087 


<0.8 


< 0.8 






Ho- 6563 


322.03 ±0.15 


306.38 ± 0.22 




248.54 ± 2.63 


[N II] 6584 


< 1.6 








[S II] 6716 


< 1.6 








[S II] 6731 


< 1.6 








[Ar V] 7006 


< 1.6 








He I 7065 


< 1.6 








[Arlll] 7136 


< 1.6 








[O II] 7320 


< 1.6 








[O II] 7330 


< 1.6 








[O II] 7333 


< 1.6 








?? 7408 


0.82 ± 0.24 








[CI IV] 7530 


< 1.6 








He II 7600 


< 1.6 








O V7611 


< 1.6 








OIV 7713 


< 1.6 








C IV 7726 


< 1.6 








[Arlll] 7751 


< 1.6 








[CI IV] 8046 


< 1.6 








CIII 8196 


< 1.6 








He II 8237 


1.76 ±0.16 








H P16 8502 


< 1.6 








H P15 8545 


< 1.6 








HP14 8598 


< 1.6 








HP13 8665 


< 1.6 








HP12 8750 


1.29 ±0.26 








HP10 9014 


1.7 : 








[S III] 9069 


< 1.6 









In particular, He n A 1640, Nm] All 50, and Cm] ,11909 are not 
seen. Table[3]gives upper limits for the latter line intensities, with 
respect to He n A 1640, as seen in the FUV spectrum. 
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Fig. 8. HST FUV nebular spectrum spectrum of TS 01, not cor- 
rected for reddening, showing the Hen A 1640 A, Nv A 1240 A, 
Niv]il488A, and Civ ,11550 A lines. The fluxes are in 10~ 14 
erg cirT 2 s _1 A _1 . 

Table 3. Observed HST UV line fluxes, relative to He II 1640 



ion lambda 


bright lobe 


faint lobe 


combined 


He a A 1640 A 


1.56" 


1.15" 


(2.7 +/- 0.3)" b 


Nv,U240A 


0.31 


0.51 


0.47 +/- 0.04 


O iv /1 1402 A 


<0.06 




< 0.06 


N iv] ,11488 A 


0.17 


noisy 


0.12+/- 0.05 


Civ,U550A 


1.10 


1.33 


1.28 +/-0.10 


Nm] ,11750 A 


<0.1 




<0.1 


Cm] ,11909 A 


<0.1 




<0.1 


" The flux of He H/11640 is in units of 10 14 erg cm 2 s 1 
b All the flux in He n ,11640, including both lobes 



2.3.2. FUSE data 

The observations of TS01 with the Far Ultraviolet 
Spectroscopic Explorer ( FUSE) and their reduc tions (done 
by G.T) are described in iTovmassian et alj d2004l) . No emis- 
sion lines were detected in the observed wavelength region 
between 900 and 1200 A, except the H Ly/3 line. For the 
photoionization modelling of the nebula, it is important to 
determine upper limits to the intensities of nebular lines 
expected in this wavelength range. We proceeded in the fol- 
lowing way. From a previous model of TS01 we took the 
computed nebular continuum. We superimposed on it the lines 
Cm ,1977.020, Nm ,1989.799 and Hen ,1992.4 with FWHM of 
0. 1 A (which corresponds to the measured expansion velocity 
30 km sec -1 , see Sect. 13.21 1. To this we added the central star 
model mentioned in Sect. 12.2.11 The resulting spectrum was 
processed through the interstellar hydrogen absorption simulator 
(http://violet.pha.jhu.edu/~gak/fwebsim.html) to be compared 
with the observations. It turned out that the corresponding lines 
start to be detectable in the resulting spectrum when the line 
flux reaches approximately 7.5 x 10~ 14 ergs cirT 2 s -1 . Indeed, 
the wavelength region that we are exploiting here is very 
complicated. Apart from the different interstellar absorptions 
and terrestrial airglow, the lines in this region also lie at the 
edges of the detectors where they overlap and errors are much 
higher compared to other regions to the red. 
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Fig. 9. LH infrared spectrum of TS01 between 22 and 28 /mi, 
showing the [O iv] /125.89 /im and [Ne v] /124.32 /im lines. 



2.4. Mid-infrared spectroscopy 
TS01 was observed using the 



iHouck et al.l 


2004 


( Werner et al. 


2004) 



on board the Spitzer Space Telescope 



observations used the Short-High (SH: 9.9-19.6 pm; R-600) 
and Long-High (LH: 18.7-37.2 /mi; R~600) modules. The aper- 
ture of the SH module is 4.7" x 11.3" and of the LH one is 
11.1" x 22.3", so the entire nebular flux was measured. The de- 
tails of the performed observations are shown in Table Q] For 
LH we used four exposure cycles of 240 s each for on-source 
and off-source observations, while for SH only on-source ob- 
servations were performed with six exposure cycles of 480 s 
each. The starting points for our interactive data reduction were 
the co-added 2D flat-fielded BCD (Basic Calibrated DATA) im- 
ages (one for each node position; pipeline version 15.3 for SH 
and 17.2 for LH). The rogue pixels were removed using the 
IRSCLEAN toofl with the aggressive parameter equal to 0. 
Then the data were processed (full extraction, trimming, defring- 
ing and averaging over c ycles) into a single spectrum per node 
position using SMARlQ dHigdon et al.ll2004 . A similar proce- 
dure has been applied for LH off-source observations and the 
obtained spectra have been subtracted from the on-source data 
for the corresponding node position, to cancel out the sky back- 
ground. The resulting spectrum between 22 and 28 pm is shown 
in Fig. [9] For the high-resolution SH module no background sub- 
traction was done since no sky measurements were taken and the 
SH slit is too small for on-slit background subtraction. Finally, 
the spectra obtained for both modules were averaged over two 
node positions and the detected nebular lines measured within 
SMART. 

The resulting intensities are listed in Table [4] together with 
the estimated uncertainties. These uncertainties do not include 
possible calibration errors. It is generally considered that the ab- 
solute flux calibration has an accuracy of 20-30%. This will be 
taken into account in the modelling (see Sect. |4.21 i. 



3 This tool is available from the Spitzer Science Center website: 
http://ssc.spitzer.caltech.edu 

4 SMART was developed by the IRS Team at Cornell University and 
is available through the Spitzer Science Center at Caltech. 



ion lambda 


flux 


uncertainty 


rad. velocity 




10" 21 W cm" 2 


10" 21 W cm" 2 


km s -1 


[Oiv]/l25.89/im 


0.38 


0.03 


-180 


[Nev] A24.32/im 


1.50 


0.20 


-140 


[Nev] A 14.32 fim 


0.82 


0.10 





Table|4]also lists the blue-shifts of the lines. One can see that 
they are consistent with the optical measurements of Tovmassian 
et al. (2001). 



2.5. XMMdata 

TS 01 has been also observed in the X-r ays, with XMM. The dat a 
acquisition and analysis is presented in Tovmassian etTajl d2009l) . 



3. What else do we know about TS01 and its 
exciting star? 

3.1. Extinction 

TS01 suffers onl y little extinction. U sing the observed Hy/HB 
and H5/HB ratios. lRicher et al.l d2002l) had found E(B - V) ~ 0.3 
mag. However, this estimate was made without considering the 
underlying stellar absorption in the Balmer lines. Due account 
for this effect significantly reduces the estimated E(B - V), as 
noted by Tovmassian et al. (2004). The extinction can also be 
estimated by considering the spectral energy distribution of the 
stellar core as observed in the far UV by FUSE. Assuming a tem- 
perature of 120 kK for the central star, Tovmassian et al. (2004) 
obtained a good fit to these observations for E(B - V) = 0.045 
mag, when using non-canonical value for Ry of 2 .3 and the in- 
terstellar reddening tables from iFitzpatrickl (Il999l) . Such a low 
value of Ry, as compared to the standard one of 3.1, was con- 
sidered compatible with the location of TS01 well outside the 
galactic disk, since the intervening dust is likely composed of 
smaller grains than in the spiral arms. However, we now know 
that the temperature of the star do minating the UV conti nuum 
is much cooler (see Sect. 13.31 and [Tovmassian et alll2009l) . im- 
plying that a steep reddening law is no t needed, after all . In the 
remainder of the paper, as w ell as in iTovmassian et all (120091) 
we use the IFitzpatrickl (fl999) reddening law parametrized with 
Rv = 3.1, and take E(B - V) = 0.03 mag, which satisfactorily 
account for the observed H and He line ratios as well as the ob- 
served continuum. Note that the absence of an absiorption dip at 
2200 A imposes an upper limit of 0.06 for E{B - V). 



3.2. Expansion velocity 

The expans i on ve locity of TS01 has been measured by 
iRicher et al.l d2003f) . This parameter is useful to estimate the ex- 
pansion cooling in the nebula. It also allows one to have an idea 
of the nebular dynamical age. We adopt v exp = 30 km s _1 . 



3.3. The stellar core 

Our understanding of the stellar core of TS01 has evolved 
considerably since the first paper where it was suggested 
that SBS 1 150+599 A is a high excitation planetary nebula 
(Tovmassian et al. 2001) with a central star having an effec- 
tive temperature above 100,000 K. Spectroscopic variations in 
the course of one single night, reported in Tovmassian et al. 
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(2004), indicated the presence of a double system with a com- 
pact star. Photometric obs ervations then unambig uously deter- 
mined a period of 3.92 h dNapiwotzki et alj|2005[) . Analysis of 
the light curve indicated that the visible star is likely an elon- 
gated ellipsoid irradiated by a source of higher energy. It also 
supported the previous conclusion that the companion must be 
a (pre -?) white dwarf or a neutron s tar. Finally, X-ray observa- 
tions dTovmassian et al.l 120071 120081) obtained with the XMM- 
Newton satellite revealed directly the light from the companion, 
which turns out to be a hot compact star! Thus, as will be shown 
later, the "cool" star is the one visible in the optical and the UV 
and it provides most of the ionizing photons. But it is the "hot" 
star which gives rise to the high excitation lines observed in the 
nebula. The best fit to the total spectral energy distribution of 
the binary core indicates that the cool component has a tem- 
perature T c ~ 60 kK while the hot component should have 7^ 
~ 170 kK. However, the determination of the temperature of the 
hot compone nt is not very accurate. Note that, in the scenario 
developed bvlTovmassi an et al (120091) . the hot component is an 
old white dwarf, which has a 170 kK temperature not because it 
is still early on its cooling path, but because it was heated by nu- 
clear burni ng of the accreted materi al on its surface. For the cool 
component lTovmassian et all d2009l) obtains the following : T c = 
58000 K + 3000 K, log g c ~ 5.1. The lower limit on the tem- 
perature is the intrinsic temperature of the star, the upper limit 
corresponds to the zone heated by irradiation. It is important to 
note that the cool component is not spherical and has not only an 
inhomogeneous temperature distribution on its surface but also 
an uneve n gravitational accelerati on. Its total luminosity is esti- 
mated bv lTovmassian et al d2009l) to be L c = 1 700 L with about 
30% uncertainty. In the following, for the sake of simplicity, we 
will consider that the cool star is sufficiently well represented by 
a stellar model atmosphere with T c = 55 kK and log g c ~ 5.1, 
with a total luminosity of 1700 L . 

The abundance analysis performed by T.R. on the cool star 
gives 12 + log He/H = 10.95 and 12 + log C/H = 7.20, with an 
uncertainty of about 0.3 dex, and upper limits 12 + log N/H < 
6.92 and 12 + log O/H < 6.80. 



3.4. Interstellar absorption of nebular UV lines 
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Fig. 10. Stellar model atmosphere fitting of carbon lines in the 
stellar core of TS 01. Observation is in black. 
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Fig. 11. Stellar model atmosphere fitting of nitrogen lines in the 
stellar core of TS 01 . Observation is in black. 



In the course of his stellar atmosphere analysis, T.R. noted 
that the observed C iv /1 1550 A and N v A 1240 A absorption lines 
were stronger than predicted by his best models. He suggested 
that these lines are probably affected by interstellar absorption. 
In that case, the intensities of the C iv A 1 550 A and N v A 1 240 A 
nebular lines are also affected by absorption. Since these lines 
are crucial for the determination of the nebular abundances in 
TS 01, we here explain how we corrected for this effect. 

We use the following notations (all quantities are a function 
of wavelength): Fs: flux extracted at the position of the star; F^: 
flux extracted at the adjacent position in the nebula; F„: real stel- 
lar flux; F ne t,: real nebular flux; F s k y : sky emission and nebular 
continuum. The optical depth due to interstellar absorption is de- 
noted r. 

We have: 



F s = (F» + F neb + F sky ) exp(-r) 



and 



F N = (^neb + F sky ) exp(--r), 
so that, in the spectrum analyzed by R.T., we have: 

F s - F N = F, exp(-r). 

Concerning the C iv /11550 line, reading out from Fig.[10l we 
find F„ exp(-r) = 8.7 x 10 -15 (black line in the figure), leading 
to exp(-r) = 0.75. Therefore, if the measured nebular flux is 
Fn eb exp(-r) = 1.28 x 2.7 x 10~ 14 = 3.45 x 10~ 14 (last column 
of Table [3]), the nebular flux after correction for absorption is 
F neh = 4.65 x 10~ 14 with an uncertainty of about 20%. 

Concerning the Nv/11240A line, from Fig. QT| we find: 
F, = 2.1 xl0~ 14 (red line in the figure), F t exp(-T) = 1.9xl0~ 14 
(black line in the figure), leading to exp(-r) = 0.9. Therefore, if 
the measured nebular flux is F ne b exp(-r) = 0.47 x 2.7 X 10~ 14 = 
1.27 x 10" 14 (last column of Table[3]), the nebular flux after cor- 
rection for absorption is F ne b = 1.41 x 10~ 14 . Here, the uncer- 
tainty is larger, since the line is weaker. We adopt 30%. 

4. Photoionization modelling 

4.1. Global strategy 

The chemical composition of TS01 can only be determined 
through photoionization modelling, since we have no direct elec- 
tron temperature diagnostic. With the observational data now at 
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Fig. 12. Density structure of the nebula. Left: The chosen den- 
sity structure of the model along the polar axis and along an axis 
perpendicular to it; Middle: The resulting Ha surface brightness 
distribution along the same axes (continuous lines), compared 
with the observed distribution (dotted lines); Right: The theoret- 
ical Ha image. 



hand, we are able to confine the range of possible abundances 
much better than in previous studies. In this paper, we try to 
make the best use of all the observational constraints. The first 
aspect concerns the morphology. The HST image (see Fig.[T|) has 
an elliptical shape, with two distinct narrow lobes. The fact that 
the Ha surface brightness distribution shows a hole in the centre 
indicates that those lobes are real (possibly due to jets) and not a 
thin disk seen edge on. Because the presence of these lobes can 
affect the ionization structure of the nebula, we choose to carry 
out the photoionization modelling with a code that allow s us to 
deal with such geometries: Cloud y _3D dMorissetll2006l) . based 
on Ferland's ID code CLOUDY dFerland et al.lll998l) . We use 
version c07.02.01 of CLOUDY and version 594 of Cloudy _3D. 

We assume that the nebula is axisymmetric with its large axis 
in the plane of the sky, and that the lobes have a circular cross- 
section. By trial and error, we choose a density structure so as to 
reproduce the observed Ha surface brightness distribution. The 
chosen density law along the polar axis and along an axis per- 
pendicular to it is shown in the left panel of Fig. [12] The resulting 
Ha surface brightness distribution along the same axes is shown 
in the middle panel of this figure (continuous lines), and com- 
pared with the observed distributiorQ (dotted lines). The right 
panel shows the theoretical Ha image, which can be compared 
with the observed image shown in Fig. [1] especially as regards 
the width of the polar lobes. Note that the density contrast be- 
tween the lobes and the main body of the nebula is very modest: 
only a factor ~ 2. The density distribution is parametrized by no, 
the value of the density at the center. For each run, we choose 
«o so that, within a circle of radius 2.5 ", our model returns an 
H/3 flux of 2.5 x 10 -14 erg cirT 2 s , whi ch corresponds to th e 
observed extinction-corrected value from iRicher et alj d2002lfl 
The value of «o is thus dependent on the distance d for which 
the computations are made. 

The distance d, in turn, results from a fitting of the theoretical 
optical/UV continuum to the observed one (taking into account 
nebular continuum, aperture effects and reddening). 

For the stars, we use model atmospheres computed by T.R. 
with the Tubingen NLTE Model Atmosphere package (TMAP). 
For the cool component, we use models tailored for our ob- 
ject. For the hot component, in absence of sufficient observa- 
tional constraints, we chose among the complete flux tables for 
H-Ni models with halo composition (May 2001) downloaded 
from http :// astro. uni-tue bingen. de/ ~rauch Those models are de- 
scribed in Rauch (2003). 



3 For comparison with the model, we have symmetrised the observed 
nebular surface brightness. 
6 Note that the models are not ionization bounded. 



4.2. Putting ultraviolet, optical and infrared fluxes on the 
same scale 

After a model has been run, the extinction-corrected line inten- 
sities are computed for each of the observing slits and compared 
to the observations. This is the best way to deal with aperture 
corrections, in particular when combining UV and optical, or IR 
and optical data. Indeed, such a procedure accounts for the ion- 
ization structure of the object under study. 

Absolute calibration of spectroscopic observations is notori- 
ously difficult. We intercalibrate the UV/optical data by forcing 
the measured value of the He n /il640/He n /14686 ratio to the 
one predicted by our photoionization models in the correspond- 
ing slits. The value of /(STIS), representing the factor by which 
the measured UV fluxes have to be multiplied in order for the 
He ii /11640/He n ,14686 ratio to be in agreement with the model, 
lies between 0.90 and 0.92 in our models. The value of /(STIS) 
is larger for models with larger electron temperature. To allow 
an easier comparison between models and observations, we fix 
the value of /(STIS) to 0.91. 

For Spitzer-IRS observations, we multiply the observed 
fluxes by a factor /(IRS) which adjusts the observed values of 
[Ne v] ,124.3 yt/m/[Ne v] /13426 (after reddening correction) to the 
one predicted by the photoionization model in the correspond- 
ing slit. The values of /(IRS) range between 0.87 and 0.95 for 
the models we considered. It might be judged unreasonable to 
scale infrared fluxes using the [Ne v] A24.3 yum/[Ne v] ,13426 ra- 
tio. However, in the electron temperature domain relevant for 
TS 01, this ratio does not vary very strongly (from Tq— 20 kK to 
40 kK, it decreases by only a factor 2). In any case, this is the 
only option we have to link together the Spitzer line fluxes with 
the optical ones, since our Spitzer data contain no H or He lines. 
Of course, in the discussion, we bear in mind this difficulty. To 
remove the model dependance of the IR fluxes correction, we fix 
the value of /(IRS) to 0.91. 

The fact that both /(STIS) and /(IRS) are found very close 
to unity is remarkable and means that the flux calibration of the 
STIS and IRS LH spectra of TS 01 is excellent. 

4.3. Judging a model 

In order to judge a model, it is convenient to divide the line ratios 
to be fitted in different categories: 

- Ratios of hydrogen lines or of helium lines: they probe the 
reddening law, the stellar underlying absorption, and the re- 
combination line theory. 

- Ratios of two different ions of the same element, such 
as [Oiv]i25.9//m/[0 in] ,15007, N v il240/[N iv] ,11486, 
[Ne v] i3426/[Ne in] ,13869, and [Ne v] ^3426/[Ne iv] 
,14720. They basically test whether the ionization structure 
is well reproduced by the model. In this category, we add 
the He n ^4686/H/? line ratio, which is more dependent on 
the ionization level of the nebula than on the abundance of 
helium. 

- Ratios of lines used to determine the chemical composi- 
tion: [Oiv]A25.9fim/Hf3, C iv ,U549/HjS, N v ^1240/H/3, 
and [Ne v] A3426/H/3. We also consider [O m] A5007/H/3 
(although it is redundant with [Oiv] A25.9pm/H/3 once the 
ionization structure is reproduced). 

Note that, in the case of TS 01, the only ratio for direct plasma 
diagnostics (i.e. electron temperature and density) that is avail- 
able is [Nev] A24.3 p.m/[Ne. v] ,13426. Unfortunately, the two 
lines come from measurements in different apertures and with 
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Fig. 13. Graphic chart to compare the reference model with the 
measured line ratios of TS01. The numbers on the top of the 
panels indicate the values of k{0) for the line ratios listed at the 
bottom. The line ratios presented in this figure were computed 
for the following pairs of slits (using the same nomenclature as 
in Table g): (0/0), (4/4), (3/0), (3/4), (1/1), (0/0), and (0/0) for 
the ratios testing the degree of ionization, and (0/0), (3/4), (1/0), 
(1/0), and (0/0) for the ratios used to determine the abundances. 



different observing techniques, and, as mentioned above, there 
is a priori some uncertainty in the relative calibration of the two 
wavelength domains. However, the fact that we find /(IRS) close 
to unity argues that the electron temperature of the [Ne v] emit- 
ting zone in our models is not far from the true one. 

For all the observables considered (usually line ratios), we 
compute the value of 



k(O) = (logO mo d - logO obs )/r(<9), 



(2) 



where O mo d is the value returned by the model, <9 b s is the ob- 
served value, and t(O) the accepted tolerance in dex for this 
observable. For each observable, the value of t(O) is chosen a 
priori considering the observational error bar, including the un- 
certainty due to reddening, and the expected ability of our model 
to reproduce a given observable. The value of t(O) is defined as 
follows: 



t(O) = log(l + AO/O), 



(3) 



where AO is the absolute value of the maximum "acceptable" 
error on the observable. We then judge our models by looking at 
their outputs presented in graphical form (see an example in Fig. 
[T3I ). A model is fully satisfying only if each of the values of k(O) 
is found between -1 and +1, and if the computed line intensities 
satisfy the upper limits for undetected lines. Of course, a prelim- 
inary condition for a model to be considered is that it returns the 
correct value of the H/3 flux, as explained in Sect. 14.11 



4.4. The reference model 

Here we present our reference model, R. This is the model for 
which all the values of k(0) are as close as possible to zero, 
taking the following characteristics for the cool star: r c =55kK, 
log g c =5.1, L c = 1700L o , and applying the extinction correc- 
tions with Ry = 3.1 and E(B - V) = 0.03 mag as explained in 
Sect. 13.11 The model contains graphite grains (as expected for 
a carbon rich planetary nebula, which is the case of TS01 as 
seen below). The grains have a standard size distribution and a 
total dust-to-gas mass ratio of one tenth the standard value. A 
larger abundance of grains would bring the predicted continuum 
around 24 /mi into conflict with the observation. In the follow- 
ing, we explore deviations from the reference model which still 
account for the observational data. 

The reference model has 7h=170kK, log gh-6.7, and a total 
luminosity Lh =2564 Ljl «o=181 crrT 3 and the following abun- 
dances, in units of 12 + log (X/H): He=10.95, C=7.84, N=7.15, 
0=6.82, Ne=6.83, S=5.65, Ar=4.70. 

The models and observations to which they are compared are 
presented in Table [5] Column 1 of this table lists the line identi- 
fications, Column 2 characterizes the observation using the fol- 
lowing nomenclature: for CFHT 2003, 1 for STIS, 2 for FUSE, 
3 for Spitzer, 4 for SDSS, 5 for CFHT 2001 . For lines which be- 
long to a wavelength range that was not observed the number 
6 is attributed. Column 3 lists the observed reddening-corrected 
line intensities (or their upper limits), in units of H/3=100 in the 
corresponding apertures. Column 4 lists the acceptable relative 
error AO/O used to compute k(O). In the case of HST, FUSE and 
Spitzer data, we estimate the value of Hfi in the relevant aperture, 
basing on our models (since, as explained in Sect. [4] they deliver 
a smoothed version of the observed surface brightness distribu- 
tion). The top rows of column 5 of the table list the character- 
istics of the reference model. The predicted line intensities in 
the relevant aperture are given in the following rows, in units of 
H/?=100 in the same aperture. For easier analysis, the next rows 
list a few important line ratios, the intensity of each line being 
measured through the aperture corresponding to the observation. 
In order to shorten the table, we do not list the lines for which the 
predictions from all our models give values smaller than 0.001 
of Hy6 (we note that such is the case for all the recombination 
lines of elements C, N, O). We list only the strongest H and He 
lines (we checked that the weaker H i and He n lines always give 
\k(0)\ < 1 - 1.5 in our models, implying that the correction for 
stellar absorption and reddening is satisfactory). 

The graphical representation of the line ratios predicted by 
model R and used to estimate the chemical composition of TS 01 
is shown in Fig. [13] This is the kind of chart that was used in 
practise when judging the models that were run. A "best model" 
is one for which all the diamonds fall as close as possible to ordi- 
nate 0. In any case, an acceptable model should have all line ra- 
tios represented by a diamond between the two horizontal lines, 
which represent a one-sigma deviation from the observed value. 
In addition, acceptable models should not return line intensities 
above the upper limits allowed by the observations. 

Figure[l4]shows the monochromatic images of the reference 
model in various emission lines. It reveals a few interesting fea- 
tures of the model: some lines, such as [O m] /15007 and [Ne m] 
/13869 arise mainly in the lobes. Other lines, such as C iv /11549, 
N v A 1240, and [Ne v] /13426 line come from the entire neb- 
ula (in agreement with what Fig. [3] suggests for [Ne v] /13426), 



7 The mass of the hot stellar component in the models depends on 
the value assumed for the gravity, which is not well constrained. For the 
photoionization modelling, what really matters is T h and L h . 
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Fig. 14. Monochromatic images of the reference model in var- 
ious lines (the values of the wavelengths are in A if they are 
larger than 900, and in /im otherwise). The x and y values are the 
coordinates in pixel units of the models. 



while O vi /11032 and [Ne vi] A7.6 //m (the latter not in the ob- 
served wavelength range) come from the innermost regions. We 
can also see that the two C m lines, C m] ,11909 and C m A977, 
although produced by the same ion, come from different regions: 
C in ,1977 has an important component coming from the central 
main body (see Fig. PHI , where the very high electron temper- 
ature allows for its excitation even if C ++ is not very abundant 
there. 

Figure [15] compares the observed energy distribution from 
TS01 with the computed one. The top panel shows the 
reddening-corrected flux computed for model R, in erg cirT 2 s _1 
A -1 , in the wavelength range 900A-40/im. The observations 
are superimposed in various colours, as indicated in the caption. 
One can see that the model reproduces the observed spectral en- 
ergy distribution quite well (except for the IRS SH observations 
which could not be corrected for sky emission, as explained in 
Sect. |2.4t . The bottom panel shows the energy distribution in the 
soft X-ray range: the blue curve is the hot star, the green triangles 
are the XMM observations. The flux from the star has been cor- 
rected for the nebular absorption (computed by CLOUDY) and 
for the interstellar absorption, taking a hydrogen column den- 
sity of 1.6 x 10 20 cirT 2 . For each computed model we checked 
that the ionizing flux does not violate the observed stellar emis- 
sion up to 200 e V. At higher energies, the observed emission may 
have another origin than the stars we consider, but it does not af- 
fect our model fitting, since we have no relevant observational 
constraints (the ion with highest ionization potential observed is 
Ne 4+ , which has an ionization potential of 97. 1 eV). 

Figure [l6]compares the energy distributions of the two stars 
considered in the modelling: the "cool" star is in red, the "hot" 
one is in blue. The sum of the two is in black. As mentioned 
in Sect. 13.31 the cool star dominates in the optical range (a few 
eV), and until about 20 eV, but it is the hot star which provides 
the photons with energies above 40 - 50 eV. So that it is the 
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Fig. 15. Comparison of the reddened spectrum of the reference 
model with observations. Top: From the UV to the IR. The 
reddened model is in black. The colour code for the observa- 
tions is as follows. Magenta: FUSE; blue: HST; green: SDSS ; 
red: Spitzer (the SH observations could not be sky-corrected). 
Bottom: X-ray domain. The model (with extinction applied) is 
in blue. The XMM observations are represented by triangles. 



cool star which provides most of the H ionizing photons, but 
it is the hot star which provides the photons responsible for the 
presence of the He n, N v, [O iv], [Ne iv], [Ne v] and [Ar v] lines. 
This is a very uncommon situation, perhaps a unique case among 
planetary nebulae: TS01 has two ionizing stars! This explains 
why our previous attempts to model the object were facing the 
difficulty that the nebula needed plenty of photons of energies 
above 54.4 eV, while the Balmer absorption lines in the stellar 
continuum indicated a moderate temperature. 

From Figs.QjJEI and Table|5] one can see that our reference 
model fits all the observational constraints very well. The only 
exception is that of the C m] A 1909, which is slightly above the 
upper limit we gave to the STIS observation. However, we con- 
sider this result to be still acceptable, since upper limit fluxes for 
unobserved lines are difficult to estimate accurately. The abun- 
dances of C, N, O, and Ne in the reference model are, respec- 
tively, 1/3.5, 1/4.2, 1/70, and 1/11, and for S and Ar < 1/30 of 
the Solar values given bv lAsplundet~ai1 ( l2005l) . 



5. The chemical composition of TS01 

5. 1 . Range of abundances for the reference model 

We now investigate the error bars on abundances that are due 
only to the uncertainties in the observed line intensities. Since 
C, N, O, and Ne contribute negligibly to the energy budget, 
it is straightforward to estimate from model R the minimum 
and maximum abundances corresponding to the minimum and 
maximum values of [Oiv] ,125.9 /mi/H/3, C iv ,11549/H/?, N v 
/H240/H/3, and [Ne v] /13426/H/3 without changing the ionization 
structure of the nebula. However, one has also to consider the 
error bars on line ratios that constrain the ionization structure: 
[Orv],l25.9//m/[0 m] ,15007, N v ,11240/[N iv] ,11486, [Ne v] 
^3426/[Ne m] ,13869, and [Ne v] ^3426/[Ne iv] ,14720. The min- 
imum values of the C, N, O, and Ne abundances in TS01 are 
obtained by a model with the lowest ionization compatible with 
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Fig. 16. Spectral energy distribution (normalized to arbitrary 
units) of the radiation from the ionizing stars in the reference 
model. Red: the "cold" star; blue: the "hot" star; black: the sum 
of the two. 

the observations and the lower limits of [O iv] /125.9 fim/Hj3, C iv 
,11549/H/?, N v ,ll240/H/3, and [Ne v] /13426/H/?. Such a model, 
Mi, is reported in column 6 of Table [5] It is derived from the 
reference model R by lowering the values of Lh, and decreas- 
ing the values of the C, N, O, and Ne abundances. With similar 
considerations, one can construct a model Ma, which will give 
the maximum C, N, O, and Ne abundances. This model, with a 
higher Lh and same Th is listed in column 7 of table [5] 

The resulting limits on the abundances of C, N, O, and Ne in 
the gaseous phas^l of TS 01 are thus: 

7.64 < 12 + log C/H< 8.05 
7.00 < 12 + log N/H < 7.32 
6.63 < 12 + logO/H<7.13 
6.76 < 12 + log Ne/H < 6.90 

The limits on the C/O, N/O and Ne/O ratios are obtained by 
considering tailored models reproducing the extreme values of 
the observed intensities of C, N and Ne lines. They are: 

0.83 < log C/O < 1.21 
0.12 < log N/O < 0.54 
-0.16 < log Ne/O < 0.17 

To derive the lower limit on He/H, one must consider the 
model with the lowest ionization compatible with the observed 
[Oiv] ,125.9/mi/fO in] A50Q1, N v ,U240/[N iv] ,11486, [Ne v] 
,t3426/[Ne m] ,13869, and [Ne v] ,13426/[Ne iv] ,14720, and the 
lowest He n /14686/H/3. The upper limit on He/H is obtained 
with similar arguments. The corresponding models, HeMi and 
HeMa, respectively, are listed in columns 8 and 9 of Table 
The resulting limits for He/H obtained in this way are 0.095 and 
0.081. In other words, the precision on the He/H abundance is 
not very good, despite the fact that we have been fitting the He n 
/14686/H/? ratio within 4% (the formal uncertainty in this ratio 
is 2% for the 2003 CFHT data, but comparison with SDSS data 
led us to adopt a higher value for the tolerance). With a more ac- 
curate upper limit on He i /15876/H/3 we could reduce the error 
bar on the helium abundance (as a matter of fact, model HeMi 
slightly violates the present upper limit on He i /15876/H/3). But 



8 The contribution of grains to the abundance of carbon is discussed 
inSect. U23l 



the uncertainty in He/H will remain larger than the uncertainty 
in the He n /14686/Hy6 ratio, mainly because the electron temper- 
ature gradient in this nebula is large and the ratio of emissivities 
of He ii /14686 and H{3 slightly varies with temperature. 

5.2. Additional sources of abundance uncertainties 

In this section, we discuss how reasonable variations of the pa- 
rameters that were so far fixed in the modelling procedure af- 
fect the derived abundances. We also discuss some more general 
problems that may have an influence on the estimated chemical 
composition of TS 01. To save space, the models that were con- 
structed to discuss these additional uncertainties are not listed 
in the paper. We will only mention their impact on the derived 
abundances. Note that all of those additional models have the re- 
quired angular size and total H/3 flux, and their abundances have 
been chosen so as to fit the observed emission line ratios. 

5.2.1 . The effect of changing the description of the stars 

So far, we have kept the parameters of the cool star fixed. Even 
if they are rather well de termined, as explained in Sect. I3.3l and 
iTovmassian et all J2009), it is important to see the effect that 
a change in those parameters implies on the derived chemical 
composition of the nebula. It turns out that an increase of 5 kK 
in T c induces a decrease in the C, N, O and Ne abundances by 
0.05-0.08 dex. A change in g c by 0.1 dex leaves the abundances 
of the fitted nebular model unchanged. 

We then explored the effects of changing the parameters of 
the hot star. By increasing the temperature of the hot star by 
lOkK, (implying a slight decrease of its luminosity in order to 
fit the observations) one decreases the C, N, O, Ne abundances 
of the fitted nebular model by about 0.2 dex. 

We also explored the effect of changing the model atmo- 
sphere of the hot star. One extreme case is to consider a model 
atmosphere composed only of H and He, instead of the Galactic 
halo chemical composition. Consequently, there are no absorp- 
tion edges in the atmosphere above 54 eV. The model which fits 
the observations presents an intense [Ne vi] ,17.6 |im emission 
(unfortunately outside the wavelength range of our IRS observa- 
tions of TS 01). Its Ne abundance is higher by 0. 15 dex than that 
of the reference model, while the abundances of the remaining 
elements are almost unchanged. 

5.2.2. Dust issues 

Concerning extinction and reddening issues, changing E(B - V) 
and Ry within limits compatible with the observed Balmer 
decrement and the observed stellar energy distribution does not 
change the abundances derived for TS 01 significantly. 

Our reference model has a dust-to-gas mass ratio of 10 _1 
times the canonical value, with the canonical grain size distri- 
bution as stated in Sect. 14.41 The chemical composition of the 
grains -pure graphite- is dictated by the fact that the object is 
undoubtedly carbon-rich, as seen in Sects. [331 and |5. II The total 
abundance of grains in the reference model is chosen so that the 
predicted infrared flux arising from the heated grains does not 
exceed the observed IRS LH flux and that it produces no signifi- 
cant dip at 2200 A since this is not observed. The total amount of 
carbon locked in grains in the reference model is 0.4 times that 
of the abundance of carbon in the gas phase. This means that the 
total abundance of carbon in the nebula (gas plus grains) is larger 
by about 0.15 dex than given in Sect. 15. II 
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5.2.3. The role of morphology 

While constructing our reference model (and all the models de- 
scribed before), we have chosen a geometry that reproduces the 
observed Ha surface brightness, including the lobes. It is inter- 
esting to experiment a simpler model without any lobes, in which 
the averaged surface brightness is the same as in the reference 
model. The abundances in such a model differ insignificantly 
from those of the reference model. We have to confess that we 
were somewhat surprised by this result, since as shown in Fig. 
IT4l the emission in such lines as [O m] A50Q1 or [Ne m] /13869 
traces the lobes very distinctly. On the other hand, one has to 
remember that the density contrast between the lobes and the 
ambient medium is only a factor 2, as seen in Fig. IT21 

5.3. Caveats 

5.3.1. The problem of Ha 

One of the intriguing problems in the observations of TS01 is 
the behaviour of the Ha line. As seen in Table [2] the observed 
Ha/H/3 ratio varies among data taken during different runs and 
at different telescopes. Since the ratios of all the remaining hy- 
drogen lines look normal, within the error bars, we are inclined 
to think that this Ha problem has no influence on the derived 
chemical composition. Nevertheless, we feel it important to try 
to understand the reason for the observed values of Ha/H/3. 

In the present study, we have done the computations with 
the full treatment of hydrogen as offered by CLOUDY (this, 
and not case B, is actually the default option in CLOUDY). 
Under the physical conditions in this nebula, one indeed does 
not expect the Balmer lines to be emitted under case B, not 
even with the added effect of collisional excitation. The ioniza- 
tion parameter of the emitting regions is high and the nebula 
is opticall y thin, making it a good candidate for case C as de- 
scribed bv lBaker & Menzefl dl938i) and reconsidered bv lFerlandl 
d!999t) . In such a case, absorption of Lyman photons from the 
star contributes to the emission of the Balmer lines, and the 
Balmer decrement depends on the number of respective Lyman 
line photons in the star. However, we are far from reproducing 
the Ha/H/3 ratios observed in the various slits. Of course, the 
computed Balmer decrement strongly depends on the fluxes at 
the wavelengths of the H Lyman lines in the model atmosphere 
used. But the differences in the Ha/H/3 ratios in the different ob- 
serving runs make it doubtful that simple stellar fluorescence can 
explain the observations. 

The reference model predicts a ratio of about 2.81. The dif- 
ferences in the observed Ha/H/3 ratios cannot have a nebular ori- 
gin since the associated time scales are far too long. 

Water vapour absorption near Ha is far too weak to explain 
the variations^. Now that the nature of the binary central star is 
better known, we can also discard the possibility that much of the 
Ha emission comes from an accretion disk. Active mass trans- 
fer in the system has ceased and, even if there is a stellar wind 
or weak remnant of accretion disk around the hot component, it 
cannot have a big influence on emission lines, since we detect 
fairly symmetric underlying absorption lines from the cool com- 
ponent at all orbital phases. These symmetric lines also imply 
that extra emission from the irradiated face of the cool compo- 
nent does not contribute significant Ha emission. 

T he remaining option is atmospheric refraction (iFilippenkol 
119821) . since the slit was not oriented at the parallactic angle for 



9 http://www.astrossp.unam.mx/sitio/abs_telurica_english.htm 



many (though not all) of the spectra with Ha/H/3 ratios differ- 
ing significantly from 2.81 (rather, usually east-west). What is 
odd, a priori, if atmospheric refraction is responsible, is that the 
lines from H5 to H/3 are observed with constant intensity ratios. 
Simulations in which we convolve the quantum efficiency of the 
slit camera used at SPIvfl with the object's very blue spectrum 
indicate that the effective wavelength is between 4000A and 
4500A. Thus, the effective wavelength, which is what is used 
to centre the object in the slit, is between the blue lines, so atmo- 
spheric refraction has very little effect upon them. As a result, 
Ha should be the only optical line that may be significantly af- 
fected by atmospheric refraction. Also, compared to the usual 
assumptions, the wavelength baseline over which atmospheric 
refraction operates is unusually large in this case, of order of 
2000 A or more. Tests using the SPM4 dataset dRicher et al.l 
12002b . in which this issue can be studied in greatest detail, 
clearly implicate the effect of atmospheric refraction since the 
spectral shape of the central star's continuum varies as a function 
of the difference between the slit position angle and the parallac- 
tic angle. Therefore, we are inclined to attribute the variations 
observed in the Ha/H/3 ratio to atmospheric refraction. 

5.3.2. Atomic data 

As noted by iPequignot & Tsarnisl d2005l) . the atomic data on 
which photoionization models are built are not of perfect ac- 
curacy. All the models we have computed rely on CLOUDY 
c07.02.01. It is not excluded that future advances in atomic 
physics, especially in the calculation of recombination coeffi- 
cients for highly ionized species, might affect the computed ion- 
ization structure. However, the fact that we now have observa- 
tional data (or stringent upper limits) on several ions of each of 
the elements C, N, O and Ne in TS 01 makes us confident in the 
robustness of the chemical composition that we have derived. 
The relatively large error bars we are obtaining on the abun- 
dances (principally due to the lack of a direct measure of the 
electron temperature in the nebula) imply that the uncertainties 
in atomic data, including the collision strengths of the lines used 
for the diagnostics, should be negligible in the total error budget. 

5.3.3. Dynamical effects 

ISchonberner et al.l d2005l) have drawn attention to the possible 
importance of dynamical effects in the thermal balance of neb- 
ulae. They make the point that the role of dynamical expan- 
sion in the cooling budget increases as the metallicity decreases. 
We have therefore included the effect of expansion cooling in 
CLOUDY, by introducing a wind cooling contribution in the rou- 
tine cool_eval . cpp: 

dynamic s.dDensityDT = (float) (2 . *fudge(Q)) ; 
CoolHeavy.expans = 

dense . pden*phycon. te*B0LTZMANN* dynamics . dDensityDT ; 

with the user defined parameter "fudge" being related to the ex- 
pansion velocity and the outer radius of the nebula by fudge= 

^exp/^out- 

All the models presented above have been computed with 
an expansion velocity of 30 km s , corresponding to the ob- 
served value (see Sect. 13. 2\ . We have tried other values for v exp 
in the equation above, but noted no significant changes in the 



10 http://www.astrossp.unam.mx//Instruments/bchivens/camrend/manual- 
english.pdf 
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output between and 200 km s _1 , the extre me values we tried. 
This r esult is at variance with the finding by ISchonberner et~aT] 
d2005l) that expansion cooling significantly reduces the tempera- 
ture with respect to a fully static model of same density structure. 

In our models, the dominant cooling process is collisional 
excitation of H Lya, and, at the ionization level predicted by 
the model, it is clear that expansion cooling must be negligible, 
unless the velocity of the jet is of the order of 1000 km s _1 . 

Could it be that the lower temperature found by 
ISchonberner et al.l d2005l) in fully dynamical models with re- 
spect to hydrostatic ones, which they attribute to expansion cool- 
ing, is actually the result of some other process? The only idea 
that comes to mind is departure from ionization equilibrium. 
For an average temperature of 30 kK, and an average density 
of about 200 cm -3 , the recombination time for hydrogen is of 
about 10 3 yr. From the apparent size, expansion velocity and 
distance to TS01, one can estimate an expansion time of ~ 
7x 10 3 yr. Therefore, the nebula should not be far from ionization 
equilibrium. On the othe r hand, the dynamical model shown in 
ISchonberner et al.l d2005l) was for a 0.595M Q star with an effec- 
tive temperature of lOOkK, corresponding to an evolution time 
of ~ 5 x 10 3 yr. The average density of the nebula in their sim- 
ulation is then about 100 ctrT 3 . In such a situation, the nebula 
is farther from ionization equilibrium. Since their star is in a 
phase where the number of ionizing photons increases with time, 
the ionization level of the dynamical model should be smaller 
than that of the corresponding static model. Hence, Lyman alpha 
cooling should be more important, and the electron temperature 
lower than in the hydrostatic model, which is indeed what their 
dynamical model gives. In TS01 the dynamical effects on the 
ionization and temperature of the nebula should be much smaller 
than in the case computed by Schonberner et al., if noticeable at 
all. In their model, the temperature drop due to dynamical effects 
is of about lOkK. Given the argumentation above, we consider 
that any dynamical effect on the electron temperature in TS01 
would be of 2-3 kK at most, with respect to the temperature we 
compute in our model. As an experiment, we computed a model 
where we use the CLOUDY parameter cextra with a value of 
2Q-20.3 gr g cm 3 s -i to simulate an extra cooling factor that re- 
duces the average electron temperature by about 3 kK with re- 
spect to the reference model. We then adjust the abundances so 
as to reproduce the observed line ratios. We find that the abun- 
dances in this model are not very different from those of model 
R. In particular, the abundance of O is not changed. The reason is 
that the model must reproduce the [Ne v] /13426/[Ne m] /13869 
ratio, which is nearly temperature-independent, and that ratios 
used to constrain the oxygen abundance ([O iv] A25.9fj.mfH/3 and 
[O m] /15007/H/3) are not very sensitive to the temperature above 
30 kK. 

5.4. Wrapping up 

In summary, considering all the possible sources of uncertain- 
ties, and adding in quadrature the various independent errors, 
we find that the elemental abundances in the gas phase of TS 01 
are as listed in Table [6] 

An additional amount of carbon, about 40% of the total el- 
emental abundance, is locked up in dust grains. Allowance for 
this component raise the carbon abundance in the nebula to 12 + 
log C/H = 8.00+0.3. 

The error bars on the derived abundances may seem large, 
when compared to the typical error bars in other PNe. However, 
one must remember that the analysis of TS 01 is much more dif- 



ficult, due to the absence of direct temperature diagnostics and 
to the weakness of the lines from metals. 

The abundances derived for the nebula are consistent with 
those derived by T. R. for the atmosphere of the cool star, see 
Sect. 13. 31 except for carbon whose abundance is larger by 0.8 dex 
in the nebula. This agreement is remarkable, given the difficulty 
of the analysis. Whether the discrepancy between the carbon 
stellar and nebular abundances is real should be examined in 
more detail. 

6. Understanding the observed chemical 
composition 

6.1. TS 01 compared to other Galactic halo PNe 

Figure [T71 shows the abundance pattern of TS01 with respect 
to those of other PNe located in the Galactic halo as derived 
by iHoward et al.1 (1 19971) . The panels in the top rows of this fig- 
ure show the elemental abundances with respect to oxygen, as a 
function of O/H, a common way to display abundances in PNe. 
As has been noted before, PNe belonging to the Galactic halo 
display a large dispersion in their abundances relative to oxy- 
geio In this respect, TS01 is not an exception. But it is ex- 
treme in its value of O/H which is significantly lower than in 
other PNe of the Galactic halo. On the other hand, the C/O, N/O, 
Ne/O ratios are similar to the highest ones found in those ob- 
jects. In the panels of the bottom row, the abundance ratios are 
computed with respect to argon and displayed as a function of 12 
+ log Ar/H. The reason for doing this is that the abundance of Ar 
is not expected to be modified with respect to the initial chemi- 
cal composition out of which the progenitor star was formed. Of 
course, the determination of the argon abundance is much less 
accurate than that of oxygen. In the case of TS 01, we even have 
only an upper limit. Yet, the bottom row of Fig. [l7]confirms the 
impression that, regarding the chemical composition, TS01 is 
different from the remaining halo PNe especially because of its 
remarkable low oxygen abundance. It may have the lowest Ar/H 
as well, but we do not know, since we have only an upper limit 
on its abundance. 

As can be seen from Fig. \T7\ the progenitors of many of the 
PNe belonging to the Galactic halo underwent considerable nu- 
clear processing, which affected not only the helium, carbon and 
nitrogen abundances in the nebulae, but also the oxygen one. 
In such a situation, the metallicity of the progenitor -generally 
identified with the oxygen abundance in the nebula- cannot be 
easily determined. If one uses the sum of the mass-fraction abun- 
dances of C, N, O, Ne, S, and Ar as a proxy for the upper limit of 
the metallicity, we find that, in Solar units, TS 01 has a metallic- 
ity of at most 1/12 of Solar, well below the upper limits for the 
remaining halo PNe. If we take the argon abundance as a proxy 

11 The dispersion in S/O, however, is likely significantly affected by 
important errors in the sulfur abundances. 



Table 6. Nebular abundances of TS 01, in various units 





12+log X/H 


uncertainty 


X/H 


mass fraction 


He 


10.95 


+ 0.04 


8.91 x IQ- 1 


2.63 x 10"' 


C 


7.84 


±0.30 


6.92 x 10" 5 


6.11 x 10" 4 


N 


7.15 


+ 0.25 


1.41 x 10" 5 


1.46 x 10" 4 


O 


6.82 


±0.33 


6.61 x 10" 6 


7.79 x 10" 5 


Ne 


6.83 


±0.30 


6.76 x 10" 6 


9.96 x 10" 5 


S 


< 5.5 




<3.16x 10" 7 


< 7.45 x 10" 6 


Ar 


<4.5 




<3.16x 10" 8 


< 8.38 x 10" 7 
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10.9 11 11.1 
12 + log He/H 

Fig. 17. Comparison of TS 01 (filled red circle) with other Galactic PNe in abundance ratios diagrams. The different black symbols 
represent: K648 (circle), DdDm-1 (square), PN006-41.9 (triangle), NGC 2242 (diamond), NGC 4361 (plus sign), M2-29 (asterisk), 
H4-1 (star), PN 243.8-37.1 (filled triangle), and BoBn-1, which actually belongs to the Sagittarius dwarf spheroidal galaxy (x sign). 
The Sun is represented by the Solar symbol. 



for the metallicity, we find that the metallicity of TS 01 is less 
than 1/30 of Solar. The two other halo PNe with very low metal- 
licity, using this criterion, are BoBn-1 (1/20 of Solar) and K 648, 
which holds the record from this point of view (1/300 of Solar). 

6.2. What we know of the progenitor ofTS01 

Let us first summarize the main features of a possible evolution 
of the stellar core of TS01 as explained in Tovmassian et al. 
(2007, 2008, and 2009). The total mass of the system is close to 
1.4 M . The mass of the cool component is (0.5-0.6) M . Then 
the mass of the hot component must be (0.8-0.9) M . For 
0.001 this implies that its progenitor had a mass of < 3 M 
dHurlev et al.l 120001) . Star formation in t h e halo ceased at least 
lOGyr ago, see e.g. iMarfn-Franch et al.l d2009l) . Since the life- 
time of a (2 - 3) M star with Z 0.001 is < 300 Myr, the first 
Roche lobe overflow occurred at least 10 Gyr ago. It resulted in a 
common envelope with inefficie nt accretion, hardly more than ~ 
0.01 M (lYungelson et al.ll995l) . The initially more massive star 
of the system turned into the currently hot component. Since we 
observe a PN now, the companion of the hot component must 
have completed its evolution in 10 Gyr. The least massive star 
able to do this had a mass of about 0.9 M after the first com- 
mon envelope stage. When this low mass star evolved off the 
main sequence, a symbiotic system formed. Mass transfer dur- 
ing this stage was negligible as well. The symbiotic stage ended 
when the low mass star filled its Roche lobe on the AGB. Again, 
common envelope formed, reducing the separation of the com- 
ponents and diminishing the orbital period to 3.9 hr. The nebula 
we observe now is the remnant of this second common envelope. 

6.3. Can standard single AGB models explain the observed 
abundance pattern? 

We consulted the available yields from recent standard evo- 
lution models for singl e s tars at low metall icities (see 
iKarakas & Lattanzioll2003l and lCristallo et al.ll20Q9h . The corre- 
sponding nucleosynthetic predictions for low-mass AGBs show 
the signature of recurring third dredge-up episodes that increase 
the surface abundances of C as well as (but to a lower extent) 
that of O and N. This is at odds with the abundance patterns de- 
termined in TS01, which exhibits an extreme O depletion. As 



a matter of fact, only massive AGB models (Mj n j > 4M ) pre- 
dict O depletion as a result of hot bottom burning (HBB) dur- 
ing the thermal pulse AGB phase (TP-AGB). Besides being in 
conflict with the binary analysis recalled above, these massive 
AGB models are facing other difficulties. First, they predict that 
carbon also should be depleted by the CNO-cycle in the enve- 
lope. Second, such massive AGB stars are predicted to expe- 
rience deep second dredge-up that increases the He abundance 
well above the observed value (number ratio of He/H ~ 0.093, 
thus Y ~ 0.26 in mass fraction). 

6.4. Towards an explanation of the chemical composition of 
TS01 

We note however that those models use standard assumptions. 
In particular they do not account for rotation-induced mixing 
that is known to affect stellar evolution and nucleosynthesis (see 
iDecressin et al.ll2009bl and references therein). Additionally, the 
rotational transport itself may have been streng thened by the 
stella r coupling i n the binary system (see e.g. Ide Mink et al.l 
l2009t iPolsetall 1 19971) . Unfortunately binary stellar models 
treating the effect of mixing on the nucleosynthesis of low-mass 
stars are not yet available, but we can use single rotating AGB 
models to infer the impact of mixing during the pre-AGB phase 
on the surface chemical composition. 

We thus compare the abundance obtained for our nebula to 
standard and rotating mode ls of a IMq star computed wi th the 
code STAREVOL (V2.92) dSiess et al.lkoOOT: ISiessll2006l) from 
the pre-main se quence to the end of t he TP-AGB phase (for 
more details, see IDecressin et al. I l2009ah . In the rotating model, 
an initial rotation of 100 km s _1 is assumed on the zero age 
main sequence. This is the mean value given bv lGaigd(ll993h 
for a star with the effective temperature of our models near the 
ZAMS (~7000 K). This observational value is obtained from 
stars in the Hyades, which are more metal-rich than our model. 
If low-metallicity stars are born with the same amount of angu- 
lar momentum we could expect, due to their compactness, an 
even higher initial velocity. Thus our model indicates the mini- 
mal chemical changes we would obtained when rotation is taken 
into account. 

The transport of angular momentum and chemical s pecies 
is driven by meridional circulation and shear turbulence dZahnl 
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Fig. 18. Comparison of the chemical composition in TS 01 with 
the results of stellar evolution models for a star with initial mass 
1M , and initial composition as indicated by the black dot. The 
model without rotation is represented in continuous line, the 
model with rotation is represented with dotted lines. The modi- 
fications of the chemical composition on the main sequence are 
indicated in green, those during 2nd dredge-up are in blue, and 
those corresponding to the TP-AGB phase in red. 



119921: iMaeder & Zahril Il998l) . The metallicity is set to Z = 
5 x 1CT 4 (i.e., [Fe/H] —1.6), which corresponds to the up- 
per limit for the metallicity derived from the Ar abundance in 
TS01. The composition is scale d to the Solar one according 
to the iGrevesse & Sauval (1 19981) mixture and enhancement in 
a-elements ([a/Fe] = +0.3 dex) is accounted for. We use the 
OPAL opacity tables dlglesias & Rogerslfl996l) above a temper- 
ature of 8 kK that acc ount for C and O enrichments, and the 
IFerguson et al.1 d2005l) data at lower temperatures. We follow 
the evolution of 53 chemical speci es from 'H to 37 C1. We use 
the NACRE nuclea r reaction rates (lAngulo et al.l 1 199*91) by de- 
fault and those by Caughlan & Fowled (Il988l) otherwise (see 
ISiess & Arnouldl2008 ). The treatment of convection is based on 
the classical mixing length formalism with q-mlt = 1-75, and 
diffusive o vershoot is add ed below the convective envelope ac- 
cor ding to Herwig d2000l) . The mass loss rate is computed with 
the iReimersi d 19751) formu la (with ?7r - 0.5) for the pre-AGB 
phase and then we use the IVassiliadis & Woodl d 19931) one dur- 
ing the AGB p hase. We also add diffusive overshoot following 
iHerwigj d2000i) to allow third-dredge up during the AGB phase. 

When comparing the results from rotating AGB models with 
those of non rotating ones, we find considerable differences in 
the predictions of the chemical composition as illustrated by Fig 
18. In particular rotation-induced mixing efficiently transports 
chemical species from the H-burning shell to the surface dur- 
ing central He-burning. This tends to decrease the surface abun- 
dances of C and O in favour of that of N through the CNO- 
cycle while He increases only mildly. As can be seen in Fig.[T8l 
rotation-induced mixing leads to low O and high N abundances 
already at the end of the second From Fig. [18] we see that in- 



troducing rotation (and thus strong internal mixing which can 
mimick binary interactions) dramatically improves the agree- 
ment of the predicted abundances with the chemical composi- 
tion of TS 01. However the O depletion in the model is consis- 
tent only with the upper limit allowed by the observations. The 
agreement would be improved if the star rotates faster during the 
phase of central H- and He-burning so that more efficient mix- 
ing induces stronger O depletion and N enhancement. As already 
explained above the chosen initial velocity for our models could 
be slightly too low for its metallicity. Furthermore an even more 
pronounced rise of the rotational rate can be the result of the first 
mass transfer when the massive primary star expands. This high 
velocity can also lead to a larger enhancement of Ne as required 
by the observations. In this way the full observational pattern of 
TS 01 can be explained, except for the He abundance which is 
still higher than observed. 

7. Conclusions 

The planetary nebula PN G 135.9+55.9 (SBS 1 150+599A), here 
referred to as TS 01, has been a source of many surprises. It is the 
most oxygen-poor planetary nebula known. Its stellar core con- 
sists of a close binary, whose period of approximately 3.9 hours 
is the shortest known of any planetary nebula nucleus. The bi- 
nary nucleus is the first known example of a double degenerate 
system in a planetary nebula. It is also the best candidate progen- 
itor of a supernova of type la. Here, we consider the chemical 
composition of the nebu lar shell in detail. In a companion paper, 
iTovmassian et all d2009) study the nature, history, and evolution 
of the binary system. 

We have obtained observational data in a complete wave- 
length range in order to pin down the abundances of half a dozen 
of elements in TS01. The data consist of imaging and spec- 
troscopy in the UV, optical and IR ranges. Optical data were 
obtained at different telescopes, resulting in high quality and re- 
liable observational material for this object which exhibits an 
exceptionally line-poor optical spectrum. 

We then defined a strategy to determine the abundances in 
TS 01. Since no direct electron temperature diagnostic is avail- 
able, it was necessary to apply tailored photoionization mod- 
elling. We did this with much care, taking into account all the 
observational information at our disposal, in order to minimize 
the uncertainties in the derived abundances. The photoionization 
modelling was performed using the pseudo-3D photoionization 
code Cloudy _3D based on CLOUDY (version c07.02.01). Using 
CLOUDY, we took advantage of its detailed treatment of the hy- 
drogen atom, which is crucial in this density-bounded nebula. 
Using the 3D features allowed us to take into account the appar- 
ent geometry of the nebula and the location, sizes and orienta- 
tions of the observing slits, thus minimizing the errors linked to 
aperture effects. The spectral energy distribution of the ionizing 
radiation was obtained from appropriate model atmospheres. 

One interesting feature of TS 01 is that both its stellar com- 
ponents contribute to its ionization: the "cool" one provides the 
bulk of hydrogen ionization, and the "hot" one is responsible 
for the presence of the most highly charged ions. Thus, one can 
say that TS 01 has two ionizing stars. This is the first such case 
discovered among planetary nebulae! 

The abundances of C, N, O, and Ne in the nebula are found 
to be respectively, 1/3.5, 1/4.2, 1/70, and 1/1 1 of the Solar value, 
with error bars of + 0.3 dex. For S and Ar, the abundances are 
less than 1/30 of Solar. This makes of TS 01 the planetary neb- 
ula with the most extreme composition known so far. In par- 
ticular, it is by far the most oxygen-deficient. The abundance 
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of helium is 12 + log He/H = 10.95, with an accuracy that is 
not as high as one might have desired (+0.04 dex). The nebu- 
lar abundances in TS 01 are found to be consistent with those in 
the atmosphere of the cool star, except for carbon which appears 
higher by ~ 0.8 dex than in the star. 

The observed properties of the stellar core allow us to put 
li mits on the masse s and a ge of the binary system, as explained 
in iTovmassian et al d2009t) and recalled in the present paper. In 
particular, the initial mass of the progenitor of the present nebula 
was ~ 0.9 M Q . We have computed state-of-the-art AGB mod- 
els with appropriate mass and metallicity, to compare with the 
chemical composition of TS 01. While anon-rotating model can- 
not explain the observed abundance patterns, the introduction of 
a 100 km s _1 initial rotation in the model (which may be seen as 
a way to mimick the mixing due to the presence of the compan- 
ion) greatly improves the comparison. Rotation-induced mixing 
leads to low O and high N abundances already at the end of 
the second dredge-up, whereas the C abundance increases after- 
wards during each third dredge-up. However to fully reproduce 
the observed pattern requires that the binary interactions (and 
especially the first episode of mass transfer) spin up the star to 
produce a very efficient mixing. This conjecture still needs to be 
verified by proper binary models. 

One feature for which we have no explanation so far is the 
fact that the carbon abundance in the nebula seems to be much 
higher than in its progenitor. This fact would however need to be 
confirmed with deeper observations in the optical and the UV. 

With respect to our first publication on TS 01 by Tovmassian 
et al. (2001), our understanding of the status and properties of 
this object has thus made enormous progress. While this object 
is unique when considering all its observed properties together, 
similar objects might be plentiful in the Galactic halo -although 
not necessarily in a stage exhibiting a detectable planetary neb- 
ula. 
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Table 5. Photoionization models versus observations. Line intensities are in units of H/3=100 in the corresponding aperture 



(1) 


(2)" 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 










R 


Mi 


Ma 


HeMi 


HeMa 






r k 


10 3 K 


170 


170 


170 


170 


170 






L h 


[L ] 


2564 


1618 


4064 


1618 


4064 










6.7 


6.7 


6.7 


6.7 


6.7 






"0 


cm~ 3 


181 


181 


181 


181 


181 


Sun'' 


















10.93 




He' 




10.95 


10.95 


10.95 


10.98 


10.91 


8.39 




C 




7.84 


7.64 


8.05 


7.64 


8.05 


7.78 




N ( 




7.15 


7 


7.32 


7 


7.32 


8.66 




O' 




6.82 


6.63 


7.13 


6.63 


7.13 


7.85 




Ne' 




6.83 


6.76 


6.9 


6.76 


6.9 


7.14 




S c 




5.65 


5.5 


5.83 


5.5 


5.83 


6.18 




Ai" 




4.7 


4.5 


4.92 


4.5 


4.92 






d 


[kpc] 


22.5 


22.3 


22.9 


22.3 


22.9 






Ft 




2.55 


2.52 


2.56 


2.54 


2.54 






M(R) 


nvir,i 

L iv± OJ 





0.14 


0.15 


0.14 


0.15 






/(STIS) 




0.91 


0.91 


0.91 


0.92 


0.91 






/'(IRS) 




0.91 


0.88 


0.94 


0.87 


0.95 








AO/0 












H8 





9.51 


0.15 


9.01 


9.02 


9.02 


8.96 


9.09 


H8 


4 


9.29 


0.15 


9.15 


9.19 


9.14 


9.12 


9.21 


H7 


() 


12.99 


(I.I 


14.18 


14.22 


14.18 


14.12 


14.28 


H7 


4 


13.6 


0.1 


14.39 


14.46 


14.36 


14.36 


14.47 


H6 


() 


25.66 


0.07 


23.64 


23.70 


23.64 


23.54 


23.81 


H6 


4 


25.71 


0.07 


23.87 


23.96 


23.85 


23.80 


24.03 


Hy 





47.84 


0.05 


47.40 


47.47 


47.34 


47.41 


47.41 


Hy 


4 


46.73 


0.05 


47.64 


47.71 


47.59 


47.65 


47.66 


Hiv 


5 


251.41 


0.04 


282.06 


282.48 


281.30 


283.21 


280.50 


Ho 


4 


311.03 


0.04 


279.89 


280.14 


279.23 


280.85 


278.43 


HeI5876 


5 


< 0.19 




0.16 


0.25 


0.11 


0.27 


0.10 


HeI7()65 


4 


< 1.44 




0.04 


0.06 


0.03 


0.06 


0.02 


Hell 1640 


1 


609.18 


0.1 


594.53 


572.27 


623.3 


607.55 


575.72 


HeII4686 


() 


77.88 


0.04 


76.53 


73.85 


79.67 


78.35 


73.60 


HeII4686 


4 


74.93 


0.04 


76.16 


74.00 


78.91 


78.52 


72.93 


HeII5412 


4 


5.43 


0.1 


6.19 


6.02 


6.38 


6.40 


5.88 


CIII]1909 


1 


< 61.44 




73.41 


89.01 


59.96 


93.05 


56.73 


CIII977 


2 


< 430.3 




104.58 


112.54 


96.78 


117.74 


92.22 


CIV 1549 


1 


1054.84 


0.2 


1050.89 


890.87 


1246.07 


916.75 


1200.76 


CIV4659 


() 


< 0.38 




0.14 


0.08 


0.26 


0.08 


0.26 


[NIII]57.2 


6 






1.23 


1.54 


0.96 


1.56 


0.93 


NIII]1750 


1 


< 60.84 




3.00 


4.29 


1.99 


4.49 


1.87 


NIII991 


2 


< 423.77 




7.54 


9.52 


5.91 


9.81 


5.69 


NIV]1486 


1 


73.05 


0.5 


65.89 


67.40 


61.95 


69.93 


58.82 


NV1240 


1 


339.3 


0.3 


330.38 


253.7 


425.38 


259.25 


410.86 


[OIII]4363 


() 


< 0.42 




0.11 


0.14 


0.12 


0.15 


0.11 


[OIII]5007 


() 


2.52 


0.3 


2.57 


3.24 


2.58 


3.35 


2.49 


OIV]1402 


1 


< 36.97 




7.68 


6.88 


10.60 


7.15 


10.05 


[OIV]25.9 


3 


13.72 


0.3 


13.57 


11.20 


20.35 


11.24 


20.13 


OVI1032 


2 


< 407.42 




50.48 


30.45 


105.06 


30.65 


102.45 


[NeIII]3869 


() 


0.79 


0.5 


0.82 


1.45 


0.44 


1.50 


0.43 


[NeIII]15.5 


3 


< 7.22 




0.21 


0.37 


0.11 


0.38 


0.11 


[NeIII]36.0 


3 


< 28.88 




0.02 


0.03 


0.01 


0.03 


0.01 


[NeIV]2424 


6 






25.20 


29.46 


20.51 


30.10 


19.93 


[NeV]3426 





57.67 


0.2 


58.32 


48.63 


66.92 


49.19 


65.84 


[NeV]14.3 


3 


29.62 


0.3 


43.25 


35.05 


50.69 


34.84 


50.91 


[NeV]24.3 


3 


54.17 


0.2 


54.44 


44.16 


63.76 


43.90 


64.03 


[NeVI]7.6 


6 






27.47 


15.98 


45.07 


15.84 


45.48 


[ArIV]4711 





< 0.15 




0.15 


0.15 


0.14 


0.16 


0.13 


[ArIV]4740 





< 0.15 




0.12 


0.12 


0.11 


0.13 


0.10 


[ArV]7005 


4 


< 1.44 




0.27 


0.23 


0.31 


0.23 


0.30 


[ArV]8.0 


6 






0.29 


0.22 


0.37 


0.22 


0.36 


[ArV]13.1 


3 


< 7.23 




0.45 


0.33 


0.57 


0.33 


0.57 


[SIII]9069 


4 


< 1.41 




0.24 


0.31 


0.17 


0.32 


0.16 


[SIII]18.7 


3 


< 7.22 




0.21 


0.25 


0.17 


0.25 


0.16 


[Sffl]33.5 


3 


< 36.1 




0.34 


0.41 


0.28 


0.42 


0.27 


[SIV]10.5 


3 


< 7.23 




6.65 


6.29 


6.88 


6.35 


6.74 


NV1240/NIV]1486 


1/1 


4.645 


0.7 


5.01 


3.76 


6.86 


3.71 


6.98 


[OIV]25.9/[OIII]5007 


3/0 


0.139 


0.4 


0.13 


0.087 


0.20 


0.085 


0.21 


[OrV]25.9/[Om]5007 


3/4 


0.194 


0.4 


0.22 


0.136 


0.35 


0.134 


0.36 


[NeV]3426/[NeIII]3869 


0/0 


73.021 


0.7 


71.33 


33.47 


151.80 


32.84 


154.26 


[NeV]3426/[NeIV]4720 


0/0 


99.039 


0.3 


127.42 


89.66 


181.34 


87.90 


185.72 


[NeV]24.3/ [NeV]3426 


3/0 


0.024 


0.3 


0.024 


0.023 


0.024 


0.020 


0.025 


[NeV]24.3/ [NeV]14.3 


3/3 


1.828 


0.3 


1.26 


1.26 


1.26 


1.26 


1.26 



" Labels for the observing slits: 0: CFHT 2003; 1: STIS; 2: FUSE; 3: Spitzer, 4: SDSS; 5: CFHT 2001; 



